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Abstract 

We study the effects of the coupHng to an Ohmic quantum reservoir on the 
static and dynamical properties of a family of disordered SU{2) spin models in a 
transverse magnetic field using a method of direct spin summation. The tendency 
to form a glassy phase increases with the strength of the coupling of the system to 
the environment. We study the influence of the environment on the features of the 
phase diagram of the various models as well as the stability of the possible phases. 

1 Introduction 

The coupling of quantum two-level systems (tls) to a dissipative environment has 
decisive effects on their dynamical properties. The case of dilute systems, in which 
interactions between the TLS can be neglected, has been extensively investigated in 
the literature ^ [2] and is now well understood. The physics that emerges when 
interactions between the tls may not be neglected has received much less attention. 



1 



In this paper we study the effects of a dissipative environment on the equihbrium 
and dynamical properties of quantum glassy systems. While spin-glass phases in 
which quantum fluctuations play an important role were found in a number of 
experimental systems 0111 El; the importance of interactions on the properties of 
tunneling defects in structural glasses was recently demonstrated 

Mean-field |Z| and finite dimensional ^ quantum spin-glass models have been 
studied in detail in the last few years. The interest in the effects of the coupling to 
the environment on the properties of quantum glassy systems is more recent. The 
quantum spherical p-spin glass model (HlEl^n) the SU{N) random Heisenberg 
model in the limit N ^ oo ^21? and the quantum random walk problem |E1 were 
studied in detail. In these cases the relevant degrees of freedom are continuous 
which makes analytical treatments possible. Here we analyze the more realistic case 
of quantum SU{2) spins with infinite range interactions between p-uplets of spins 
coupled to a bath of harmonic oscillators. For p = 2 this model reduces a model 
for metallic spin-glasses studied previously Ul| . For p > 3 the model exhibits richer 
behavior jl5l 1161 1171 118j . including the possibility of having first order transition 
lines. 

We use a method of direct spin summation (dss) first used in the context of 
disordered spin models by Goldschmidt and Lai In this method the disorder- 
averaged free-energy density is computed using the replica method to average over 
the random quenched interactions. Next, a Trotter decomposition is performed 
in order to express the partition function of the resulting single-site self-consistent 
problem as a sum over different contributions, each coming from a possible spin his- 
tory (t{t) = ±1. The continuous imaginary-time variable < r < Ph is discretized 
on a grid Tt = t(3h/M, t = 0, ■ ■ ■ , M — 1, and the partition function is computed by 
numerically performing the exact sum over the 2*^ possible discrete spin histories, 
at = CF{Tt) = ±1. Physical results are obtained repeating this procedure for various 
values of M and extrapolating to M — > oo. 

We find that the coupling to the enviromement favors the appearance of the 
spin-glass phase reducing the strength of the quantum fluctuations that tend to 
destabilize it. For p = 2 the phase transition is always second order. For p > 3 
there exists a tricritial temperature T* below which quantum fluctuations drive the 
transition flrst order. T* decreases with the strength of the coupling to the bath. 
For p > 3 a dynamic transition precedes the equilibrium phase transition. The 
coupling to the bath also stabilizes the dynamic glassy phase. 
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The organization of the paper is as follows. In Sect. [21 we introduce the coupled 
p-spin-bath model and the formalism that we use to solve it. In Sect. 01 we discuss 
the numerical method and present our results. Sect. contains our conclusions. 



2 The model 

We are interested in disordered SU{2) spin models in a transverse field described 
by Hamiltonians of the type 

TV 

Hs = HL-Tj2at (1) 

i=l 

with 

N 

Hl = - E Jiu■■■,^,K■■■K■ 

ii<-<ip 

Here, a^,a'^,a^ are the standard Pauli matrices, Jii---ip denotes a quenched random 
exchange between p spins, T is the transverse field introducing quantum fluctuations 
and is the total number of spins. The sum runs over all possible p-uplets of spins. 
The model is then fully-connected and mean-field in character. It is completely 
determined once the order p and the precise distribution V[J] of random interactions 
are chosen. We consider the case in which the random independent variables Jii^...^ip 
are Gaussian with zero mean and variance plj'^ /IN^^^ . The scaling of the variance 
with N is chosen so as to ensure a good thermodynamic limit. 

We study the thermodynamics and some aspects of the non-equilibrium dynam- 
ics of the quantum spin model coupled to a quantum environment assumed to be 
in thermal equilibrium. We model the coupling to the environment by assuming 
that each spin in the system is coupled to its own set of N/N independent har- 
monic oscillators with the total number of them. The bosonic Hamiltonian for 
the isolated reservoir is 




The coordinates xi and the momenta pi satisfy canonical commutation relations. 
For simplicity we consider a bilinear coupling, 

TV TV 
i=l 1=1 
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that involves only the oscillator coordinates. The Hamiltonian for the coupled 
system is then given by 

H = Hs + H^ + Hsb + Ha (5) 
where we added a counter-term, 




whose effect is to eliminate a possible mass-normalization induced by the coupling 
to the bath |2j. 

The partition function of the combined system for a particular realization of the 
bonds 

Z = TY[e-^^] (7) 

involves a sum over all states of the system and of the bath. The trace over the 
variables of the bath can be performed explicitly using standard techniques |21 1213 
I2T] . The result of tracing out these variables can be expressed in terms of the 
spectral function of the bath, 

N 

W ^ lj2^^Siu;-iVi)=6,,Iiu;). (8) 

We chose to study the effect of an Ohmic bath parameterized as 

I{uj) = 2r]hujQ{0Jmax " w) (9) 

with r] the friction constant, tOmax an ultra-violet cut-off and Q(x) the Heaviside 
theta- function. 

This problem can be mapped onto a classical Ising spin system using the Totter- 
Suzuki formalism |15[ 1191 1^ . This amounts to writing the path-integral for the 
partition function as a sum over spin and oscillator variables evaluated on a discrete 
imaginary-time grid, Tt = f3fi/Mt, labeled by the index t = 0, . . . , M — 1. Periodic 
boundary conditions on the discrete time-axis are imposed due to the fact that 
the evaluation of the partition function involves a trace. To recover the correct 
representation of the trace the limit M — > oo should be ultimately taken. The finite 
M expression yields a sequence of M-approximants to the asymptotic M — > oo 
formula. 
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The M-th approximant of the "reduced" partition function obtained after inte- 
grating out the bath reads 



Z = Trj^tj exp 



M-l N 



a M-L iV M-l N 

t=0 n<---<ip t=0 4=1 

M-l N 

- E E(i-^k*')C(i-o 

t,t'=0 i=l 



(10) 



where 



C 



A 
B 

{t-t') 



1 

2 
1 
2 

2r] 



sinh 
coth 



pr 

M 



cosh 



(11) 
(12) 



_., , cosh (cj/?7i((t-t')/M- 1/2)) sinh2(a;/?n/2M) , , 



a; sinh {ujj3h/2) 

The trace represents the sum over ah 2^^*^ distinct classical Ising spin config- 
urations, a\ = ±1, for each spin, i = 1,...,A^, evaluated at each time-slice 
t = 0, ...,M-1. 

The disordered averaged free-energy is calculated using the replica trick [23] 



I3F = -In Z 



lim 



- 1 



(14) 



n^O n 

After some standard manipulations, and up to some irrelevant factors, we obtain 

n M-l 



Z" = f[ Yl /"pg"*^*'pA"*''*'exp(-iVP(A,g)) 

a,b=lt,t'=0 



(15) 



with 



n M-l 



P{A,Q)= E E 

a,b=l t,t'=0 \ 



I ^1 



a2 t2 

/^atbt' '-^ //^atbt'\mp , /i xo-b i -i /^atbt' \ 

I ^ "4M2^^ )^ + C(t_t,)() (1-g , 



InTri^atie^-^^/ , 

n M-l ^atW 



H, 



n M-l 



e//= E E + ^ ^ + 

a,fe=lt,t'=0 " 



a=l t=0 



(16) 
(17) 



where the bullet is used to distinguish the ordinary power from the matrix power. 
In the thermodynamic limit, N oo, the integrals in can be evaluated with the 
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saddle point method at the expense of exchanging the — > oo and n — > hmits. 
The disordered-averaged free-energy per spin is then 



Pf 



hm 

n— >o n 



where Aq and Qo are such that, 
6P{Q,A) 



6Q 







5P{Q,A) 



?o,Ao 



6A 



. 



(18) 



(19) 



Hereafter we omit subscripts in the saddle-point values Qo and Aq. The disorder- 
averaged entropy per spin is easily obtained from the disorder-averaged free-energy 
density and reads 



s 



a,b=l t,t'=0 



n M-l 



2 j2 



2AP 



i^qatbt'yp _ ^ab^^ _ Qatbt 



dp 



+ 



sinh(2/?r/M) 



cosh(2/3r/M) --^Q 



a(t+l)at 



a=l 



(20) 



Another physical observable of interest is the magnetic susceptibility, 

dM 



X 



dh 



(211 



h=0 



where M = N ^J2iLi {o'f) is the total disorder-averaged magnetization and h a 
longitudinal external magnetic field. In terms of Q"*^*' the susceptibility is given by 



n n M-l 



(22) 



a,b=l t,t'=0 



In Eqs. (|2n|) and (|22|) the right-hand-sides should be evaluated at the saddle-point 
values. 

The matrix elements Q*^*^*' are the order parameters of the model. 



Q 



atbt' 



1 

N 



E«-r)- 



(23) 



1=1 



Because of the translational invariance in the Trotter time direction the diagonal 
terms in the replica indices depend on the time-difference only. 



Q 



atat' 



qd{t - t') . 



(24) 
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Notice that due to the periodic boundary condition qd{t) = qd{M — t). In addition, 
as QdiO) = 1, only {t — t') = 1,2, ■ ■ ■ , int^ need to be considered. The ofF-diagonal 
elements in the replica indices, Q*^*^* with a ^ b, are t and t' independent as shown 
by Bray and Moore j24j . 

In order to determine the different phases of the model, we consider the following 
Ansdtze: 

Paramagnetic Phase 

The matrices Q and A are taken to be diagonal in replica space 

Q'^*^*' = qa{t - t')6''^ , A'""''' = Xd{t - t')5''^ . (25) 

Using Eqs. ((TB|) . (|T7)) and (|TH|) the disorder-averaged free-energy per spin can be 
expressed as 

o2 j2 M-l r,2 t2 M-1 

= - 1) E I'At - - ^ + E - lnTr^..}e^^7; (26) 

with 

, M-l M-l 

HlTf = 172 E - + E + ^^V(*+«) . (27) 

Here, qd{t — t') and Xd{t — t') are obtained self-consistently from the extremum 
condition by summing over all 2*^ spin configurations cr* = ±1: 



Trj^t} 



HPTf t t' 



Qdit-t') = ^ ^ , (28) 



X,{t-t') = LLPqP-^^t-t')+M%_,y (29) 

Equilibrium spin-glass phase 

In order to characterize this phase we use a one-step replica symmetry breaking 
(rsb) Ansatz, 

A-tbt' ^ ^Xd{t-t')-Xea)S''' + Xeae''' . (31) 
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Using Eqs. (jlbf) . (|T7|) and ((TH)) the disordered-averaged free-energy density becomes 



with 



TT^sg 
^eff 



4M 



— In / dA" (Tr 

m 



M-l 



M-1 



1 

M2 

M-l 



A 



ea , \/2Aeai 

2; 

M M 



M-l 

E-* 



(32) 



(33) 



and the integration measure 



dX 



dx 

=e 2 



(34) 



Here and in what fohows ah integrals over x go from — cxd to 00. 

As in the paramagnetic phase, the order parameters qd{t — t'), q^a, Xd{t — t') and 
Aea are determined self-consistently from the extremum conditions that involve a 
sum over all 2^'^ spin configurations, cr* = ±1: 

s9 \ m-l 



Qdit - t') 



Xd{t - 1') 



JdX (Tri^.je^^//) 

/ ri-esfl \ m— 2 / rj-esg 

/d^ (Tr{,.|e^=//) (Tr{,.}e^'=//Et^VM 
JdX (Tr{,t}e^^// 

^gr^(t-0+M2C(,_,), 

:^Qea ■ 



(35) 

(36) 

(37) 
(38) 



As it has been discussed in a number of papers on classical and quantum [T5| 
1161 1171 126j spin-glass models, two choices for the determination of the breaking- 
point m lead to different physical results. The use of the extremum condition, 
that corresponds to looking for the value of m that renders the disorder-averaged 
free-energy stationary, leads to 



m 



/-I 



m 



(p-l)gL + ln J dX (Tr{,.|e^' 



esg \^ m 



(39) 
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where 



(40) 



With this choice one describes the equihbrium properties of the model. 

This Ansatz yields the exact solution 17 to the spherical version of the p > 3 
model. When p = 2 the spherical model is solved by a simpler replica symmetric 
(rs) form. The stability of the one-step Ansatz for quantum SU (2) models can 
be tested by extending the analysis of de Almeida and Thouless |28| . When the 
lowest eigenvalue, that is also called the replicon, vanishes, the one-step RSB Ansatz 
becomes marginally stable. When the replicon is negative, this Ansatz is unstable. 

We do not expect the one-step RSB Ansatz to be stable everywhere in the phase 
diagram in the case of discrete spins. Indeed, by evaluating the replicon at the 
order parameters and break-point obtained from the extremum conditions we found 
that the one-step RSB Ansatz is unstable in the full spin-glass phase when p = 
2 [Sherrington-Kirkpatrick (sk)] model indicating the need to break the replica 
symmetry further. In the case of the p > 3 classical spin model, the one-step RSB 
Ansatz is unstable below a temperature Tg < Ts as shown by Gardner [211 in the 
classical case. Thus, the solution for the classical Ising p spin model also requires 
full RSB at very low temperatures. Tg depends on the parameter p and, as expected, 
it tends to = J when p ^ 2+ and it vanishes when p ^ oo. 

When quantum fluctuations are taken into account we thus expect to find a 
Gardner line of instability. A careful study of how this line depends on T and the 
coupling to the bath requires to solve the quantum problem at rather low tem- 
peratures. This is done in Sect. where we compute the location of the Gardner 
instability line. As seen in Fig. d the region where the one-step RSB static Ansatz 
is unstable is quite small. Outside this region the one-step RSB Ansatz is exact 
and can be studied to study the properties of the p > 3 quantum SU{2) model. 
Elsewhere, and for p = 2, we shall regard this solution as a suitable approximation 
to the correct solution. 

Dynamic spin-glass phase 

The marginality condition leads to a different equation for m. With this con- 
dition one requires that the saddle-point is only marginally stable, i.e. the matrix 
of quadratic fluctuations has a zero replicon eigenvalue (and one does not impose 
the condition of extreme on m). It has been checked by comparison to the real- 



9 



time dynamics jU], that this condition yields the freezing transition of the spherical 
quantum p-spin model with p > 3 coupled to the oscillator reservoir at the initial 
time t = jlDj . Here we use it as an indication of where such a dynamic transition 
line should be located for the SU{2) quantum spin systems. 

Adapting the calculation of de Almeida and Thouless 2H] to the quantum prob- 
lem under study we find that the replicon eigenvalue is given by 

\r = P-2Q + R (41) 

with 

P = l-kqP-h, 

R = -kq^ea'r. (42) 

The factors r, u and t are 

/ dX Tr^t e^'^ff TV,* e^^// (Et (^'/M) 
r = ^a-a^a^a'^) = ^ ^ , ^ ..,,^,,n ^ , (43) 



JdX ( Tr^te^-ff 

M2 



^ = 172E(^"*^''^"'^') (44) 



tT 

JdX (Tr^.e<nY ' fTr,*e<//(EtCT7M)yTr,.e<//(E«'^V*7M2) 



JdX Ti^te^-ff 



+ ^ \ " I ^at ht'ar br' \ 

* = Xf4 2^ ^ ^ I 



U'tt- 



zjdsg \ ™- 2 / Trdsg , 



JdX Tr^te^-// T^rate"-ff{Ett' f^V^ /M^ 



JdX i^i„te"-f^f 
Here, we have defined 

k^'^p{p-l). (46) 

and Hfj^j- is the Hamiltonian of by Eq. (j^HI)- Finally, the values of qd{t — t') and qea 
are fixed by the extremal conditions. 



dsg 



(45) 
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3 Results 



In this Section we describe the outcome of solving the equations we derived in the 
previous Section and we discuss how the coupUng to the Ohmic bath of harmonic 
oscillators modifies the behavior of the spin model. 

3.1 Numerical method 

The free-energy density and derived magnitudes depend on the parameter M that 
in practice takes finite values. Several strategies were proposed to study the limit 
M ^ CO. Usadel and Schmitz [201 noted that M should be such that f3T/M <^ 1. 
For low temperatures this criterium becomes quickly impractical since one cannot 
perform the complete sum over states for such large values of M. As an alternative, 
these authors proposed to use a Montecarlo procedure to estimate the sum over 
configurations when M is large \M\ . 

In this paper we use another method that has been previously used to study the 
isolated quantum SK 19 and the SU(2) p-spin models in a transverse field. We 
compute each magnitude using a direct spin sumation (dss) with M ranging form 
M = 8 to M = 13. The six values obtained are the input for a polynomial extrap- 
olation which gives the M — > oo limit for the considered magnitud. For almost all 
cases, we extrapolate using the (1/M)^ law [HII. Special care is taken for the study 
of the function (?d(r) that is known on the imaginary-time grid 0, • • • 

We interpolate these sets of points with splines for each value of M. As the value of 
M increases the density of points also increases and we have a better description of 
the curve. Finally we use the result of the function q^^ir) at r as the input for the 
polynomial extrapolation. Because of the additional difficulties introduced by the 
functional nature of (/^(t) we considered a linear extrapolating scheme in this case. 

Before discussing the behavior of the model in detail let us illustrate the de- 
pendence on the number of time slices M used in the DSS numerical method. In 
Fig. n we show the diagonal order parameter ^^(t) as a function of r at constant 
temperature and transverse magnetic field for the p = 3 model. Six values of M are 
shown, M = 8, • • • , 13. In Fig. [21 we display the free-energy density of the different 
phases of the p = 2> model with rj = 1. The four curves correspond to three values 
of M, M = 8, 9, 10, and the result of the extrapolation to M — > oo. In both cases 
the variation of the curves with M is indeed very smooth. 

Even if the extrapolating method of Goldschmidt and Lai IT9l is simple to im- 
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0.2 0.4 0.6 0.8 

rlfih 



Figure 1: Upper panel: The diagonal imaginary-time dependent order parameter 

qd{T) in the p = 3 model as a fimction of T/{j3K) — tjM with t the Trotter time, t = 
0, 1, • • • , M. The temperature is T = 0.3 and the transverse magnetic field is F = 0.8. 
The coupling to the bath isr] = 1. The solutions for M = 8, 9, 10, 11, 12, 13 are shown 
from top to bottom, (with lincpoints, the points being the numerical data and the 
lines representing the result of splines) and the lowest curve (thin line) is the result 
of the extrapolation to M — > oo. Lower panel: The limiting curve limM^oo 9d^(''") 
for three couplings to the environment: 77 = (bottom), 77 = 0.5 (middle) and rj = 1 
(top). The other parameters are the same as in the upper panel. (See the text for 
the details of the extrapolated method used.) 
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Figure 2: Disorder-averaged free-energy density, /, of the quantum SU{2) p — i 
model at temperature T — 0.3 as a function of the transverse magnetic field F. The 
coupling to the bath is ?7 = 1.0. The three phases of the model are represented: a 
physical paramagnet labelled PMl, an unphysical paramagnet that one discard on 
physical grounds labeled pm2 and the spin-glass. The values of / obtained for finite 
M are shown with thin hues [M = 8 (bottom), Af = 9 (middle) and M = 10 (top)]. 
The result of the extrapolation to M ^ oo is displayed with bold lines. 

plement and very efficient, it is also limited when the temperature decreases. The 
extrapolation is less clear at lower temperatures. 

3.2 Static phases 

Let us first discuss the effect of the environment on the static phase diagram of 
the p = 2 (sk) and p > 3 quantum SU{2) spin models. The spin-glass disorder- 
averaged free-energy densities are obtained using the one-step RSB Ansatz discussed 
in Sect. El We discuss the limits of validity of this Ansatz below. 

As in other disordered quantum spin models |15[ I16| \17\ I22j two paramagnetic 
solutions coexist. As in the spherical p-spin model coupled to a bath the one labeled 
pm2 in Fig. [2 can be discarded since its entropy becomes negative at sufficiently 
low temperatures. Thus, we do not discuss it further in this paper. 

The critical line {Ts,Ts) separating the paramagnetic (pm) and spin-glass (sg) 
phases is determined by the values of the pairs (T, F) where the physical paramag- 
netic (called PMl in Fig. [21) and spin-glass free-energy densities cross. 

We show in the upper and lower panels of Fig. |31 the static critical line in the 
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(T, r) plane separating a high T, high T paramagnetic phase (pm) from a low T, 
low r spin-glass (sg) for the p = 2 and p = 3 models, respectively. The three curves 
in each figure correspond to (77 = 0) and two non zero couplings (ry = 0.5, 1) from 
bottom to top. 

For both models, the classical transition temperature, Tf"'^'^, corresponding to 
Ts — > 0, remains unchanged by the coupling to the quantum heat reservoir. This 
value is Tg = J for when p = 2 and it coincides with the one given by Gross 
and Mezard, Ts ~ 0.67, for the classical problem with p = 3 [63\ . 

For the three values of r], the static critical transverse field, TsiT), is a decreasing 
function of T, which is consistent with the fact that quantum fluctuations tend to 
destroy the glassy phase. We also see from the figures that the coupling to a 
quantum thermal bath favors the formation of the glassy phase: the coupling to the 
environment effectively reduces the strength of the quantum fluctuations that tend 
to destroy it. For any value of the temperature that satisfies T < Tf"'^'^ the extent 
of the spin-glass phase is larger for stronger couplings to the bath. Moreover, we 
observe that the effect of the bath is stronger for lower temperatures. 

When p = 2 the transition is always continuous and second-order thermody- 
namically. For p = 3 instead, as in the spherical case EI E] and the isolated 
quantum SU{2) model jl5j . an interesting change from a second-order to a first- 
order transition appears. We demonstrate these statements by displaying in Figs.|3] 
and El the behavior of the free-energy density, entropy and susceptibility of the 
p = 3 SU (2) spin model as a function of the transverse field for T = 0.5 > T* and 

T = 0.3 < t;. 

At sufficiently high temperatures, T > T*, one finds a spin-glass solution for 
increasing transverse fields until the break-point m reaches the value m = 1. The 
values (T, F) where m = 1 coincides with the ones obtained by analyzing the crossing 
of the free-energy densities of the paramagnetic and spin-glass solutions. Thus, for 
the chosen temperature T > T* this is the critical transverse field. Even if the 
Edwards- Anderson parameter, qea, and the diagonal element, (/^(r), are non-zero at 
this point in parameter space, one can check, as shown in Fig.^ that the entropy and 
susceptibility do not show a jump. Thus, for T > T* the transition is discontinuous 
[due to the jump in Qea and but of second order thermodynamically. 

The situation is different at lower temperatures. In Fig. [51 we show the free en- 
ergy, entropy and susceptibility of the p = 3 SU{2) model for T = 0.3 < T* . In this 
case, the point in which the free-energy of the paramagnetic and spin-glass solution 
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Figure 3: Upper panel: Static phase diagram for the p = 2 model as obtained using 
the DSS technique for finite number of time slices M and extrapolating the data to 
M — > oo. The three lines correspond to rj = 0,0.5,1, from bottom to top. Lower 
panel: Static phase diagram for the p = 3 model obtained using the same numerical 
method. The continuous line (dashed line) indicates a second order (first order) phase 
transition. The critical lines continue below the lowest value of T for which we trust 
the algorithm, T fa 0.25, to reach a quantum critical point at T = 0. 
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Figure 4: Free energy density, /, entropy, s, and susceptibility, X) as functions of 
the transverse field, F, for the p = 3 model at T = 0.5 > T* for three values of the 
coupling to the bath, q = 0,0.5,1. The continuous (dashed) hue corresponds to the 
paramagnetic (glassy) phase. The entropy and susceptibility are continuous at the 
transition indicating a second order phase transition 
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cross corresponds to m < 1 and as shown in the figure this leads to a discontinuity 
of the entropy and susceptibility. In this case, the transition is discontinuous and 
first-order thermodynamically. 

In Fig. El we show the dependence of qea and m on the critical temperature Tg 
for three values of the coupling to the bath. The model is again the p = 3 quantum 
SU{2) spin-glass. As already mentioned we observe that for all temperatures Qea is 
different from zero, leading to a discontinuous phase transition, m equals one for 
Tg > T* but m < 1 for Tg < T*. The figure also shows that T* decreases with 
increasing coupling to the bath r/. Again, this result is reminiscent of what found 
in the spherical case |lUj . 

3.3 Stability of the one-step static solution 

In order to study the stability of the one-step solution we evaluated the replicon 
eigenvalue Xpt on the values of the order parameters and m obtained from the static 
solution, and we searched for the parameters (T, F) such that vanishes. In the 
classical limit this yields Gardner's classical critical temperature that takes a rather 
low value, Tg{T = 0) 0.25 [22]. Since we expect to find a decreasing value of 
the instability temperature with the strength of the transverse field, we need to 
control the numerical algorithm for T < 0.25. Even if this might seem, at first 
sight, impossible, we managed to obtain sensible results keeping reachable values 
of M, M < 13, since the small values of the transverse field compensate the large 
value of P in the condition (3T /M <^ 1. 

First, we analyzed the p = 2 case that corresponds to the SK model in a tranverse 
field. In the absence of the environment we found that the one-step RSB solution is 
not stable in the full spin-glass phase supporting the idea that the solution to the 
statics of this model needs a full RSB scheme, just as in its classical limit and in 
contrast to recent claims in the literature |34j . 

In Fig. [7|we compare the static critical line {Ts,Ts) as found from the one-step 
RSB Ansatz, with Gardner's line of instability for the p = 3 model. We see that the 
region where the one-step RSB static Ansatz is not stable is quite small. Since we 
only trust the extrapolation from low values of M to M — > oo above temperatures 
of the order of T « 0.1, we do not explicitly extrapolate the instability line to lower 
temperatures. Nevertheless, the existing data suggest that in the zero temperature 
limit the static critical transverse field, F^, and Gardner's critical field, Tq, do not 
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Figure 5: Free energy, /, entropy, s, and susceptibility, Xi as function of the trans- 
verse field, r, for the p = 3 model at T = 0.3 < T* for three values of the coupling to 
the bath, r] = 0, 0.5, 1. The continuous (dashed) line corresponds to the paramagnetic 
(glassy) phase. The entropy and susceptibility are discontinuous at the transition 
indicating a first order phase transition 
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Figure 6: Dependence of m and qea on the critical temperature Ts (F, 77) for the 
p = 3 model. Three values of tj are considered, rj = 0, 0.5, 1. For increasing values of 77 
the interval in which m = 1 increases and, hence, the region where a thermodynamic 
first order transition occurs decreases. 



coincide TsiTs = 0) > Tg{Tg = 0). 

3.4 The dynamic transition 

As already explained in Sect. I^lthe value of m found by setting the replicon eigen- 
value to zero leads to different equations that encode some information about the 
non-equilibrium relaxation dynamics of the system. Using this prescription we ob- 
tained, for the p > 3 models a different critical line that lies above the static 
transition. This result is similar to those found in a series of other classical |25j 
and quantum jl51 1161 [T7j problems. In Fig. |Hlwe compare the static and marginal 
critical lines for the p = W quantum SU{2) model. We chose a larger value of p 
to make the difference between the two lines easier to visualize. The glassy static 
region is smaller than the glassy region determined by the marginality condition. 
When approaching the glassy phase from any direction in parameter space, the 
dynamic transition, associated to the line of marginal stability, occurs before the 
static one. As on the critical static line, the curve determined with the marginal 
stability criterion is made of two pieces, on one of them the transition is of second- 
order (indicated with a solid line on Fig. IS)) and on the other the transition is of 
first-order (indicated with a dashed line on the same figure). The first-order nature 
of the dynamic transition is displayed by, for instance, a jump in the asymptotic 
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Figure 7: Comparison between the static critical line {Tg^Ts) and Gardner's insta- 
bility line {Tq, Tq) for the p — i model with 77 — 0. 

value of the averaged internal energy. The marginal tricritical point occurs at higher 
temperature than the static one. 

The external noise also has a strong effect on the dynamic critical line. The 
stronger the coupling to the environment (larger value of r/), the larger the spin- 
glass region in the phase diagram. This is also shown in Fig. [H] where a couple of 
curves, corresponding to r/ = and 77 = 0.5 are drawn (see the caption in the figure 
for more details). 

Finally, let us mention that there is an empirical relation between the value of 
the parameter m as found from the marginality condition and how the fluctuation- 
dissipation theorem is modified in the real-time non-equilibrium relaxation of the 
quantum model 0^]. Using this relation and interpreting then the parameter m/T 
as an effective temperature we find that the modification of the fluctuation- 
dissipation theorem, and hence Tgjj, depend on the strength of the coupling to the 
bath. 

4 Conclusions 

In this paper we studied the effect of an Ohmic quantum bath on the statics and 
dynamics of quantum disordered SU{2) spin models of mean-field type. We found 
that the coupling to the enviromemnt favors the appearance of the spin-glass phase 
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Figure 8: Comparison between the static and marginal critical lines for the p = 10 
model. The solid lines represent second order transitions and the dashed lines first 
order ones. The set of curves with Tf"'^'' « 0.6 (thin lines) are the static transition 
lines for rj = (below) and i] = 1 (above); the set of curves with T^^'^^^ « 0.82 (bold 
lines) correspond to the dynamic transition for r/ = (below) and rj — 1 (above). 



reducing the strength of the quantum fluctuations that tend to destabilize it. As in 
the case of the spherical model |1U1 lllj the phase transition is always second order 
for p = 2. For p > 3 there exists a tricritial temperature T* below which quantum 
fluctuations drive the transition first order. T* decreases with the strength of the 
coupling to the bath. For p > 3 a dynamic precedes the equilibrium phase transition. 
The coupling to the bath also stabilizes the dynamic glassy phase. 

It would be interesting to check if the same tendency to ordering appears in 
macroscopic spin models in finite dimensions. One could attempt to study this 
problem in the context of frustrated spin magnets or the much studied, numerically 
and analytically, quantum SU (2) spin chain with and without disorder. 

This problem is of interest for possible implementations of quantum computers 
where the interaction of the system with its environment needs to be controlled. 
The effect of an environment on the properties of Griffith phases has also been the 
focus of a hot debate [36! • We expect to report on these problems in the future. 
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